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ABSTRACT 



yA/ms. The aim of the project is to contribute to the characterisation of the spectral energy distribution of the coolest brown dwarfs 
discovered to date, the Y dwarfs. 

Methods. We obtained z-band far-red imaging for six Y dwarfs and a T9+Y0 binary with the OSIRIS (Optical System for Imaging 
and low Resolution Integrated Spectroscopy) instrument on the 10.4-m Gran Telescopio de Canarias (GTC). 

Results. We detect five of the seven known Y dwarfs in the z-band, infer their optical-to-infrared colours, and measure their proper 
motions. We find a higher dispersion in the z - 7 and z - H colours of YO dwarfs than in T dwarfs. This dispersion is found to be 
correlated with H - w2. The high dispersion in the optical-to-infrared colours of Y dwarfs and the possible turn-over towards bluer 
colours may be a consequence of the presence of sulfide clouds with diff'erent thicknesses, the depletion of alcalines, and/or gravity 
eff'ects. 
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1. Introduction 



The classical stellar sequence (iMorgan et al.ll943]) has now been 
extended to include very cool objects that bridge the gap be- 
tween brown dwarfs and planets, tha nks to the advent of the 
Wide Infrared Survey Explorer (WISE: IWright et al.l201Ql) satel- 
lite. T he "Y" class, originally proposed by "Kirkpatrick et al." 
(1 19991). now consists of 15 m embers (Gushing et al. 2011; 
lLiuetal.1 12011: Kirk patrick et al.ll2012i: iTinnev et al. 2012) with 
spectral types equal to, or later than, YO (tentative classification) 
and eflTective temperatures below 500 K. Two additional candi- 
dates have been reported, but their faintness currently hampers 
spectral typing (Liu et al. 2012; Luhman et al. 201 1). 

The number of late-T dwarfs with spectral types later 
than T6 has increased dramatically over the past few years 
thanks to large-scale surveys su ch as the UKIRT Infrared Dee p 
Sky Survey (UKIDSS; Lodieu et alj|2003; [Pinfield et al. 2008; 
tBumingham et al. 2009, 2010a, b), the Canada-Erance-Hawaii 
Brown Dwarf Sur vey (Delorme et al. 2008a,b; Revle e t al.ll2010l: 
lAlbert et al. 11201 l l). and WISE ( Kirkpatrick et al. 2012). The T/Y 
transition is characterised by a rapid shift of the peak of the spec- 
tral energy distribution from the near-inf rared to the mid- inf rared 
and by narrower peaks in the H band ( Gushing et al .11201 ih with 
the possible disappearance of the p otassium absorption beyond 
0.7 microns ("Leggett et al.l l2012bl) . Because of to their nature 
and low temperature, it is important to characterise their spec- 
tral energy distributions from the optical to the mid-infrared to 
understand the chemistry at play in these cool atmospheres. 



* Based on observations made with the Gran Telescopio Canarias 
(GTC), installed in the Spanish Observatorio del Roque de los 
Muchachos of the Instituto de Astrofisica de Canarias, in the island of 
La Palma. 



Our goal is to contribute to the spectral energy characterisa- 
tion of the Y dwarfs, pushing the limits of the Gran Telescopio 
de Ganarias (GTG). This letter describes the optical imag- 
ing obtained for six Y dwarfs with OSIRIS (Optical System 
for Imaging an d low Resolution Integrated Spectroscopy; 
iGepa et aDl2000l) on the 10.4-m GTG. We present the imaging 
campaign and the associated data reduction and astrometry. We 
expose our results and place them into con text with respect to th e 
spectral energy distributions of T dwarfs (iLeggett et al.ll2012bk . 



2. Far-red optical imaging 

2.1. Observations 

The OSIRIS instrument is mounted on the 10.4-m GTG op- 
erating at the Observatory del Roque de Los Muchachos (La 
Palma, Canary Islands). The detector consists of two 2048x4096 
Marconi GGD42-82 separated by an 8 arcsec gap and operates at 
optical wavelengths, from 365 to 1000 nm. The unvignetted in- 
strument field-of-view is about 7 by 7 arcmin with an unbinned 
pixel scale of 0.125 arcsec. We employed a 2x2 binning because 
it is currently the standard mode of observations. 

We imaged six of the seven Y dwarfs published by 
iGushing et all ilOUh and the T9-FY0 binary resolved by 
'Liu etal.' ('2012") using the Sloan z filter available on OSIRIS. 
Table [T] provides the original names with the first four digits of 
the right ascension and declination as provided in Gus hing et al 
(i2011.) along with their current spectral types ( Kirkpatrick et al 
l2012h . This information is supplemented with the coordinates 
measured on the GTG OSIRIS images (when the target is de- 
tected), the z-band magnitudes with their associated error bars, 
the dates of observations, and the numbers of individual im- 
ages with their respective integration times in seconds. Eigure 
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Table 1. GTC/OSIRIS astrometry, photometry (or 3cr lower limit s), observing i nformat ion, and photometry for six Y dwarfs 
and a T9+Y0 binary. The near-infrared photometry is taken from Leggett et alJ (l2012ah . except for the //-band magnitude of 
WISEJ1828-h2650 (Kirkpatrick et al. 2012). 



WISE J. . . (SpT) 


R.A. dec 


z 


Date 


ExpT 


/ 


H 


z-J 


Ha COSd 


1^6 




hh:mm:ss.sss 


mag 


dd/mm/yy 


sec 


mag 


mag 








0146+4234 (YO) 


01:46:56.576 +42:34:09.80 


24.10+0.13 


03/09/12 


50x50 


19.40+0.25 


18.71+0.24 


5.39+0.27 


-0.52+0.08 


-0.11+0.08 


0410+1502 (YO) 


04:10:22.933 +15:02:42.91 


22.66+0.09 


15/09/12 


40x50 


19.44+0.03 


20.02+0.05 


3.22+0.09 


+ 1.20+0.08 


-2.17+0.08 


1405+5534 (YOp) 




>23.85 


06/08/12 


45x50 


21.06+0.06 


21.41+0.08 


>2.79 






1738+2732 (YO)^ 


17:38:35.585 +27:32:58.28 


22.80+0.09 


17/06/12 


30x50 


20.05+0.09 


20.45+0.09 


2.75+0.13 


+0.32+0.10 


-0.39+0.10 


1828+2650 (>Y2) 




>24.46 


09/09/12 


40x50 


23.48+0.23 


22.85+0.24 


>0.98 







2056+1459 (YO) 20:56:29.028 +14:59:54.64 23.09+0.08 15/06/12 40x50 19.94+0.04 19.96+0.04 3.66+0.09 +0.89+0.10 +0.61+0.10 
1217+1626 (T9+Y0) 12:17:57.144 +16:26:35.99 21.60+0.03 15/12/12 15x60 17.83+0.02 18.18+0.05 3.77+0.04 +1.41+0.10 -1.96+0.10 



^ We also targeted this object in the / band with OSIRIS. We measured a 3cr lower limit of / > 25.3 mag, translating into i - z> 2.4 mag and 

/ - / > 5.8 mag. 



[T] shows the GTC OSIRIS finding charts (two arcmin a side) of 
all the observed objects (circled) with north up a nd east to the 
left. Ta ble [H also lists revised / and H magnitudes ('Legge tt et al.l 
l2012ah and the associated z-J colours plotted in Fig.[2l 

The observations were conducted on 15 and 17 June, 6 and 
13 August, 3, 5, 9, and 16 September, and 15 December 2012. 
All observations were conducted under average seeing of 0.8- 
1.3 arcsec, photometric or clear conditions, grey time, and air- 
mass less than 1.5. Bias and skyflats were obtained during the 
afternoon or morning of the respective nights. For each Y dwarf, 
we obtained three to five series of ten frames with 50 sec on- 
source individual integrations, except for th e coolest (WISE 
Jl 828+2650 tentatively classified as >Y2; 'Kirkpatric k et al.1 
[2QT2h for which we repeated the series four times to try to de- 
tect the object in the z-band and for WISE J1217+1626 AB, for 
which we obtained three series of five frames of 60 sec. 

2.2. Data reduction and photometry 

We reduced the OSIRIS Sloa n z-band images in a standard man- 
ner with IRAlG](Todv 1993). Eirst, we subtracted the mean bias 
and divided by the normalised averaged master skyflat to each 
individual science frame. Then, we combined each set of ten im- 
ages taken without dithering, and finally we combined those sets, 
applying the off'sets to create a master science frame. In the case 
of WISE J1217+1626AB, images were sky-subtracted before 
the combination. 

We photom etrically calibrated the z-band images with Sloan 
standard fields ( Smith et ani2QQ2h or with Sloan images of the 
same field, using point sources with photometric errors below 
0.02-0.05 mag. We measured instrumental magnitudes for 5 to 
10 stars in those fields and derived zero points ranging from 
27.976+0.037 to 28.091+0.049 mag. 

We detected four of the six Y dwarfs and the T9+Y0 binary 
(unresolved in our images) with OSIRIS in the z-band (Tabled 
Eig.[T]). We measured the photometry of the closest source to the 
nominal position of the WISE source (Table[T]), except for WISE 
J1410+1502 and WISE J1217+1626 AB, which have moved be- 
yond the five-arcsec circle drawn in Eig. [T] due to their high 
proper motion (~2.5 arcsec/yr). To measure the magnitudes, we 
performed aperture and point-spread function (PSE) photometry 
using DAOPHOT in IRAE because of the presence of objects close 



^ IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation 



to some of our targets. We chose an aperture equal to ~3 times 
the full-width-at-half-maximum and checked that the target was 
well- subtracted by our PSE analysis without residuals. We trans- 
formed the instrumental magnitudes into apparent magnitudes 
using the zero-points derived for the photometric standard fields 
observed on the same night. The final calibrated magnitudes of 
the four Y dwarfs detected on the OSIRIS images are given in 
Table [U For the remaining sources, we list the 3(T lower limits 
(Table [TJ computed from the root-mean- square of the sky at the 
nominal position of the target compared to the peak of the flux 
of three nearby (non- saturated and isolated) stars. After averag- 
ing all images for WISE J 1828+2650, we can set a lower limit 
of z > 24.46 mag, translating into z-H>\.6 mag (and z - / > 0.9 
mag) for this (tentative) Y2 dwarf. 

2.3. Astrometry 

We astrometrically calibrated the combined science frames with 
IRAF (daof ind, ccxymatc h, and ccmap rou tines) using bright 
point sources from 2MASS (ICutri et al.ll2QQ3k or SDSS DR8 in 
the case of WISE J1405p5534. We typically found 50-100 stars 
(depending on crowding) in the OSIRIS field-of-view, resulting 
in an astrometric calibration better than 0.1-0.15 arcsec. The 
GTC OSIRIS coordinates (right ascension and declination) of 
the five objects detected in the z-band are listed in columns 2 
and 3 of Table [T] We measured proper motions, consistent with 
the values quoted by Marsh et al. (2012) and Liu et al. (20^, 
thanks to the 2.0-2.5 year baseline between the GTC and WISE 
observations. The error bars of each component of the proper 
motion take into account the WISE astrometric accuracy (0.15- 
0.17 arcsec) and conservative errors on the GTC centroid of 0.15 
arcsec (Table [TJ. We find high motions as expected for objects 
lying at distances of only a few parsecs. 

3. Discussion and conclusions 

We reported in this letter the first far-red optical det ection of 
four YO dwarfs and the unresolved T9+Y0 binary from lLiu et al.l 
(2012), which are among the coolest brown dwarfs known to 
date. All four YO dwarfs with z-band detection have magnitudes 
in the range 22.8-24.1 mag, implying z-J and z - H colours 
spanning 2.7-5.2 mag and 2.3-5.4 mag, respectively (Fig. O. 
We also presented lower limit of z - / of 2.79 mag and 0.98 mag 
for WISE J1405-h5534 (TOpec?) and WISE J1828-h2650 (>T2), 
respectively. 

Figure[5] shows the z - / and z-H colours of the YO dwarfs as 
a function of spectral type (black dots with error bars and arrows 
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Fig. 1. Finding charts for the six Y dwarfs and the T9+Y0 binary imaged in Sloan z with GTC/OSIRIS. North is up and east is left. 
Images are 1x1 arcmin Targets are ordered by increasing right ascension from left to right starting from the top left comer. The open 
circles have a radius of five arcsec and mark the WISE position of the Y dwarfs, showing the motion of the object when detected 
(small red circle with two arcsec radius). 
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Fig. 2. z - J (left) and z - H (right) colours of YO dwarfs (filled dots for detection s and arrows for l ower limits) as a function 
of spectral type comp ar ed to the colours of late-T dwarfs from UKIDSS (diamonds: iBumingham et al . 201 Ob ), Sloan (aste risks; 
iGolimowski et al . 2004: Knapp et al.l l2004l : IChiu et al.|[2006D . Canada-France brown dwarf survey (squares: lAlbert et al. 11201 ll) . and 
the sample of iLeggett et al. (2012a triangles). 



for lower limits) along with the colours of T dwarfs. Here we 
adopted a spectral type of T9.5+/-0.5 for the binary T9+Y0. The 
dispersion observed for YO dwarfs is greater than for T dwarfs. 
We add ed to this plot UGP S J07 2227 .5 1 -05403 1 .2 cl assified as 
TIO bv lLucas etaP (l2010l) and'Legge tt et al.l (l2012bl) . but pro- 
posed as the T9 standard by|Cushing et an (l201lh . 

Figure [3] displays the (z - H,H - w2) colour-colour dia- 
gram for the four YO dwarfs and the T9-I-Y0 binary (filled dots) 
detected in z and for the two sources with upper limits (left- 
pointing arrows). This plot t akes into account the improved 
near-infrared photometry fromLeggett et al. (2012a), compared 
to the original values in Table 2 of Kirkpatrick et al. (2012). 
We note a trend towards bluer colours with increasing H - w2 
from late-T dwarfs to YO dwarfs, except for WISE J0146-h4234, 



which looks peculiar in all diagrams. This trend seems to start 
off' around a spectral type of ~T8 and agrees with the bluer Y-J 
colours observed with cooler temperatures for the Y class and 
the excess of ffux in the 840-940 nm region in the spectrum of 
WISE J2056-F-1459 c ompared to UGPS J0722 27.5 1-05403 1.2 
(ILeggett etal.ll2012ah . The T9-hY0 binary from lLiu et all (l2012h 
shows an intermediate position between late-T and YO dwarfs. A 
possible explanation for the relative changes in the z-band may 
be associated to the disappearance of the strong atomic potas- 
sium and sodium bands into molecules (e.g. KCI2 or NaCl) and 
possibly other alkalis as we move towards cooler temperatures. 

The spread in the z - / and z - H colours of YO dwarfs is 
-2.5-3.0 mag, much wider than the spread observed throughout 
the T dwarf sequence (Fig. [2]) and greater than the spread in the 
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Fig. 3. (z-H,H-uj2) colour-colour diagram for our targets (filled 
dots and arrows for l ower limits), late-T dwarfs (open triangles; 
lLeggettetal.ll2Q12ah . and UGPS J0722227.5 1-05403 1.2 (open 
triangle with error bars; Lucas et al.ll2Q10l) . 
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Fig. 4. (z-H,Y-H) colour-colour diagram for our targets (filled 
dots and arrows for lower limits), UGPS J0722227.5 1-054 03 1 .2 
(open triangle with the error bars; Lucas et aP l2010l) , and 
the known T9 CFB DSIR 145829.0-F101343.0 (open triangle; 
iDelorme et aD2010h . 



J - K colours of L dwarfs (iHawley et al.ll2002]) . T his spread can 
be qualitatively explained using recent models by iMorley et alJ 
(^12), which propose that the amount of sulfide clouds present 
in the atmospheres of these Y dwarfs strongly impacts their spec- 
tral energy distributions (note that this is not the case for iron or 
silicate clouds which play a role at higher tempe ratures where 
the L/T transition takes place). Indeed, Figure 5 of iMorley et al.l 
(12012") shows that the presence of sulfide clouds decreases the 
emerging flux in the 0.8-1.3 microns with lower efl'ective tem- 
peratures (starting ofl' below -900 K). Hence, the clouds af- 
fect primarily Y and / but also the z-band, which is what we 
observe. Moreover, this efl'ect seems to increase with decreas- 
ing sedimentation efficiency (the /sed parameter i n the models 
of Marley et al. (2002)). Similarly, Figure 10 of iMorlev et"aD 
(I2OI2I) demonstrates that the combination of a wide range of 
gravities and presence of sulfide clouds with a variety sedimen- 



tation efficiencies increases the spread in the / - if colours of 
brown dwarfs with temperatures below 600 K. These sulfide 
clouds may be responsible for the spread observed in the optical- 
to-infrared colours of Y dwarfs, suggesting that the atmosphere 
of WISE JO 146 may contain thicker clouds and higher gravity 
than the other YO dwarfs in our sample. Table 3 of iLeggett et al.l 
(2012a) confirms that brown dwarfs with a temperature of 400 K 
influenced by thicker clouds would be fainter in Y than those 
afl'ected by thin clouds. As pointed out above, the same efl'ect 
could occur in the z band. Moreover, we observe a correlation 
between the z - ^ and Y - H colou rs of the three YO d warfs 
(Fig. 131) common to our sample and Leggett et al.l (l2012ah . The 
bluer objects inz- H turn out to be also rather blue inY - H. If 
this trend holds, we would predict Y - H = 0.96 mag for WISE 
J0146-h4234. 

Finally, it is important to complete the sample of z-band 
measurements for Y dwarfs in the southern hemisphere and ob- 
tain far-red optical spectroscopy, which is only available for 
WISE J2056-hl459 (iLeggett et al.ll2012ah . This may show dras- 
tic changes as indicated by the z-band photometry, which can 
probe changes in the atmospheric properties responsible for the 
far-red spectral energy distribution of these brown dwarfs . Our 
photometry would also help the modellers to reproduce this 
part of the spectral energy distribution, disentangling the ef- 
fects of clouds, gravity, and pressure-broadened potassium and 
sodium bands. Thirdly, the relatively blue optical-to-infrared 
colours of YO dwarfs and the possible turnover towards even 
bluer colours with lower temperatures could have a signif- 
icant impact on the strategies currently used to identify Y 
dwarfs in large-scale surveys, such as UKIDSS, CFBDS, Pan- 
Starrs (Deacon et a l. 2 011), and lame syn optic survey telescope 
(ILSST Science CoUaboFation et al.ll2009l) . 
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